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Abstract 

We give an analytic (free of computer assistance) proof of the existence of a classical 
Lorenz attractor for an open set of parameter values of the Lorenz model in the form of 
Yudovich-Morioka-Shimizu. The proof is based on detection of a homoclinic butterfly with a 
zero saddle value and rigorous verification of one of the Shilnikov criteria for the birth of the 
Lorenz attractor; we also supply a proof for this criterion. The results are applied in order 
to give an analytic proof of the existence of a robust, pseudohyperbolic strange attractor 
(the so-called discrete Lorenz attractor) for an open set of parameter values in a 4-parameter 
family of three-dimensional Henon-like diffeomorphisms. 
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1 Introduction 

The main goal of the paper is a proof of the birth of the Lorenz attractor in the following 
model: 

A = y, 

Y = -XY + yX{l-Z)-6X^ (1) 

Z = -aZ + /3X‘^, 

where parameters a, /3, 7 , 6 and A can take arbitrary finite values. It is well-known 
that the classical Lorenz equations 


X = -a{x - y), 
y = rx — y — xz 
z = —hz + xy 


( 2 ) 


can be always brought to form ([T]) by the transformation x = X^/2, y = V 2 + Xj , z = 

(r — 1)Z -\ - with a = b, (3 = 7 = a(r — 1), <5 = 1, A = a -|- 1. We therefore 

a a[r — 1) 

call system m the extended Lorenz model. 

Not all values of parameters in system m correspond to real values of the parameters 
(r, a, b) , so the proof of the Lorenz attractor in model ([T]) may not always imply the existence 
of the Lorenz attractor in the classical Lorenz model. In particular, the values of a and A, for 
which the existence of the Lorenz attractor is proved in the present paper (see Theorem[T]), do 
not correspond to the Lorenz model; instead, for these parameters values one can transform 
system ([1]) to the system obtained from the Lorenz model by the time reversal. 
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Nevertheless, system ([T]) is interesting by itself because it is a normal form for certain 
codimension-3 bifurcations of equilibria and periodic orbits ©[5]. Thus, our result on the 
existence of the Lorenz attractor in system ([1]) automatically proves the emergence of the 
Lorenz attractor (or its discrete analogue) in a class of such bifurcations and, hence, in a 
large set of systems of various nature. In particular, in this paper, using Theorem[T]we prove 
the existence of discrete Lorenz attractors in a class of three-dimensional polynomial maps 
(3D Henon maps). 

We consider the case 7 > 0, so by scaling the time and the coordinates we can make 
7 = 1. When q:/3 > 0, with the scaling /3 = abjB one can bring system ([T]) to the form: 

w = y, 

Y = X - XY - XZ - BX^ (3) 

Z = -a{Z - X‘^). 


In [Ij it was shown numerically that system ([3]) possesses the Lorenz attractor in an open set 
of values of (a. A) for each B G (—1/3, -|-oo). In this paper we obtain an analytic proof of the 
same fact for sufficiently large B. The idea is that for B = -t-00 equations ([3]) can be solved 
explicitlj0. Therefore for large B the system can be analysed using asymptotic expansions. 

We perform such expansion for system ([T|) with 7 > 0, d > 0. By coordinates and time 
scaling it can be brought to the form: 


W = F, 

Y = X - XY - XZ - X^ (4) 

Z = -aZ + f3X^. 


Obviously, large B in system ([3]) corresponds to small /3 in Q. The following theorem is our 
main result. 


Theorem 1 For each sufficiently small [3, in the {a, X)-plane there exists a point a = 
a(i8),X = X(/3) for which system Q possesses a double homoclinic loop (o homoclinic 
butterfly) to a saddle equilibrium with zero saddle value. This point belongs to the closure 
of a domain Vla{I 3) in the {a, X)-plane for which system Q possesses an orientable Lorenz 
attractor. 

In §[2]we give a definition of the Lorenz attractor, following [6]. This definition requires 
a fulfilment of a version of a cone condition at each point of a certain absorbing domain. 

We verify this condition by using one of the Shilnikov criteria mu- A partial case of these 
criteria was proven by Robinson PEU], however it was done under additional “smoothness 
of the foliation” assumption which is not fulfilled by system Q at small fl. Therefore we 
provide a full proof of the part of Shilnikov criteria relevant to our situation, see Theorem [2j 
In § [3] we prove Theorem [1] 

In § 0] we prove (Theorem [3]) the existence of discrete Lorenz attractors in a class of 
three-dimensional Henon maps. A discrete Lorenz attractor is a generalization of the strange 
attractor in the period map of a time-periodic perturbation of an autonomous flow with a 
Lorenz attractor unEi]. It is a pseudo-hyperbolic attractor in the sense of nan]; like in 
the continuous-time Lorenz attractor, each orbit in the discrete Lorenz attractor has positive 
maximal Lyapunov exponent, and this property is robust with respect to small smooth 
perturbations of the map. 

The discrete Lorenz attractors were found numerically in several 3D Henon maps m 
dans] and in certain models of non-holonomic mechanics [lanz]- The 3D Henon maps are 

^The same holds for B = —1/3, which case will be considered in a forthcoming paper; the case of arbitrary 
B > —1/3 remains out of reach. 
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particularly important because they serve as a zeroth order approximation to the rescaled 
first-return maps near various types of homoclinic and heteroclinic tangencies [igiisiisoiEii 
[2^ 123] , Therefore, because of the robustness of the discrete Lorenz attractor, by showing 
that a 3D Henon map has such attractor we also show that the corresponding homoclinic or 
heteroclinic tangency bifurcation produces such attractor; moreover, when the corresponding 
tangency is persistent (e.g. due to the Newhouse mechanism [Ml |25l |26l [27]) one obtains 
infinitely many coexisting discrete Lorenz attractors in a generic map from the corresponding 
Newhouse domain (see [22] for examples and more explanation). Obviously, in order to 
implement this construction completely rigorously, the numerical evidence is not enough, i.e. 
an analytic proof of the existence of the discrete Lorenz attractor in various classes of Henon 
maps is needed. Our Theorem [3| in §[4|is the first example of such proof. 

We find that a certain codimension-4 bifurcation happens at certain parameter values to 
the fixed point of a 3D Henon map and show that the normal form for this bifurcation is given 
by the period map of a small time-periodic perturbation of system (|3|) at small /3. Then, 
Theorem [T] implies Theorem [3| immediately. A less degenerate bifurcation (of codimension 
3) that also occurs in the same class of maps produces, as a normal form, system ([3|) with 
B = 0. This particular case is known as the Morioka-Shimizu system |28] . While a very 
detailed numerical analysis of the dynamics of this model is available [29l [301 Ell E21 l33] 
and the existence of the Lorenz attractor in the Morioka-Shimizu system causes no doubt, a 
rigorous analytic proof of this fact is lacking (see discussion in m)- Its absence impedes the 
further development of the mathematical theory of the discrete Lorenz attractors. Obtaining 
such proof is a challenging and important problem; in particular, it would be interesting to 
know if a rigorous numerical proof of the existence of the Lorenz attractor in the Morioka- 
Shimizu system is possible, similar to that reported in [401141] for the classical Lorenz model. 

2 Shilnikov criterion of the birth of Lorenz attrac¬ 
tor 

Let us recall a definition of the Lorenz attractor. There are several approaches to it, which in¬ 
clude classical geometric models by Guckenheimer-Williams [461147] and Afraimovich-Bykov- 
Shilnikov IMIE] and their later generalisations [481149] . Here we define the Lorenz attractor 
in the Afraimovich-Bykov-Shilnikov sense. 

Consider a smooth system of differential equations which has a saddle equilibrium state 
O with one-dimensional unstable manifold. Take a cross-section H transverse to a piece of 
the stable manifold W^{0) (for system ([4]) the cross-section will be the plane z = /Sja). The 
one-dimensional manifold W^{0) consists of three orbits: O itself, and two separatrices, r_|_ 
and r_. Let both r+ and r_ intersect the cross-section at some points M_|_ and M_. Let 
Ho be the first intersection of the cross-section with W^{0), and n+ and n_ be the regions 
on the opposite sides from Hq, so H = H^ U n_ U Hq. Assume that in the cross-section there 
is a region D C Hq which is forward-invariant with respect to the Poincare map T defined 
by the orbits of the system. When a point M approaches Hq from the side of H-i. we have 
TM ^ M_|_, and when it tends to Hq from the side of n_ we have TM ^ M_ (see Fig. [T]). 
Thus, the invariance of D implies that M_|_ G D and M_ G D. Let us introduce coordinates 
{u, v) in H such that Hq has equation u = 0, n_|_ corresponds to u > 0, and n_ corresponds 
to rt < 0. We write the map T in the form 

u = f±{u,v), v = g±{u,v), 

where “-I-” corresponds to u > 0 and ” corresponds to u < 0. The map is smooth at 
u / 0 but it becomes singular as u —)• 0 because the return time to the cross-section tends 
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Figure 1: Afraimovich-Bykov-Shilnikov geometric Lorenz model. 


to infinity as the starting point approaches the stable manifold of the equilibrium O. One 
can analyze the behaviour of the orbits that pass a small neighbourhood of O by solving 
Shilnikov boundary value problem (see e.g. formulas from [36], Section 13.8). The result 
gives us the following asymptotics 


f± = u±±A±\u\'' + o{\u\''), g± = v±±B±\u\'^+ o{\u\'^), (5) 


where («+,?;+) = M+ and (u_,u_) = M_, and ly = |Ai |/7 where 7 is the positive eigenvalue 
of the linearisation matrix at O and Ai is the eigenvalue which is nearest to the imaginary 
axis among the eigenvalues with negative real part; it is assumed that Ai is real and the rest 
of the eigenvalues lies strictly farther from the imaginary axis. We assume that the saddle 
value Ai +7 is positive, i.e. u < 1. This implies that if the separatrix values A±. are non-zero 
df± , dg± 


then 


00 , 


0 as u — 0 , i.e. the map near Ho is expanding in the u-direction and 


du dv 

contracting in the ?;-directions. The main assumption of the Afraimovich-Bykov-Shilnikov 
model in [35] is that this hyperbolicty property extends from the neighbourhood of Ho to the 
entire invariant region D. For the first return map in the form ([5]) this assumption is written 
as the following set of inequalities: 


iiiiU/^r'iio • ii/:iL < (1 - ii(/o-'iio)(i - iKiL), 


( 6 ) 


where we omitted the indices ±; the notation || • ||o = sup(„ II ' II is used. These 

inequalities are an equivalent form of the so-called invariant cone conditions, which were 
later used in [371ISH] for studying statistical properties of the Lorenz attractor, in [39l jH] 
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for the numerical verification of the match between the geometrical Lorenz model and the 
Lorenz system itself, etc. 

In [ 6 ] it was shown that conditions ([ 6 ]) imply the existence of a stable invariant foliation 
of D by leaves of the form u = h{v) where h is a certain uniformly Lipshitz function 
(different leaves are given by different functions h). One can, in fact, show that the foliation 
is C^-smooth, see galls]; however, this is not important here. The map T takes each 
leaf of the foliation into a leaf of the same foliation and is uniformly contracting on each of 
them. The distance between the leaves grows with the iterations of T, i.e. the quotient map 
T obtained by the factorisation of T over the invariant foliation is an expansive map of the 
interval obtained by the factorisation of D over the foliation. This map has a discontinuity at 
u = 0: T(+0) = rt+, T(—0) = tt_ (see Fig. [T|). The fact that the quotient map is expansive 
makes an exhaustive analysis of its dynamics possible, see [441 145] . The dynamics of the 
original map T and the corresponding suspension flow can be recovered from that of T. 

The structure of the Lorenz attractor is described by the following theorem due to 
Afr aimovich-Bykov- S hilnikov. 

Theorem [G] [35] Let V he the closure of the union of all forward orbits of the flow 
which start at the set D. Under conditions above, the system has a uniquely defined, two- 
dimensional closed invariant set A <Z V {the Lorenz attractor) and a {possibly empty) one¬ 
dimensional closed invariant set E {which may intersect A but is not a subset of A) such that 
{!) the separatrices Tj- and the saddle O lie in A; 

{2) A is transitive, and saddle periodic orbits are dense in A; 

{3) A is the limit of a nested sequence of hyperbolic, transitive, compact invariant sets each 
of which is equivalent to a suspension over a finite Markov chain with positive topological 
entropy; 

{4) A is structurally unstable: arbitrarily small smooth perturbations of the system lead to 
the creation of homoclinic loops to O and to the subsequent birth and/or disappearance of 
saddle periodic orbits within A; 

(5) when S = 0, the set A is the maximal attractor in T>; 

{6) when the set E is non-empty, it is a hyperbolic set equivalent to a suspension over a finite 
Markov chain {it may have zero entropy, e.g. be a single saddle periodic orbit); 

(7) every forward orbit in V tends to A UT, as t ^ +oo; 

{8) when E / 0, the maximal attractor in V is c/(VF“(E)) = A U VF“(E); 

{9) A attracts all orbits from its neighbourhood when An E = 0. 

Let us describe one of Shilnikov criteria for the birth of the Lorenz attractor (as dehned 
in this theorem) from a pair of homoclinic loops. We restrict ourselves to the symmetric 
case; namely, we consider a C'’'-smooth (r > 2 ) system in n > 2, and assume that it is 

invariant with respect to a certain involution R. 

Let the system has a saddle equilibrium O with dimlT^(O) = n — 1 and dimiy“(0) = 1. 
We assume that O is a symmetric equilibrium, i.e. RO = O. With no loss of generality we 
assume that the involution R is linear near O |5r)j . Denote the characteristic exponents at 
O as Ai,..., An and 7 , such that 


7 > 0 > Ai > maxReAj 

i >2 

(in particular, 7 and Ai are real). The eigenvectors cq and ei corresponding to 7 and, resp., 

Ai must be ii-invariant. We assume that Rcq = —cq and Rei = ei; in particular, the two 

^In this was proved for C^-small perturbations, however the argument works for (7°“-small perturbations 
without a significant modification. A similar result is also in IHIIT]. 
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unstable separatrices are symmetric to each other, r_|_ = RT-. 


Homoclinic Butterfly. Assume that both unstable separatrices r_|_ and r_ return to O 
as t — 7 > +00 and are tangent to the leading stable direction ei (it follows from the symmetry 
that they are tangent to each other as they enter O). 

Saddle value. Assume that the saddle value a = 7 + Ai is zero (hence, the saddle index 
u = IA 1 I /7 equals to 1 ). 

As the homoclinic orbits F-i- are tangent to each other, we can take a common cross- 
section n to them near O; we can make 11 symmetric, i.e., RH = IT. As we mentioned, the 
cross-section is divided by W^{0) into two parts, 11+ and n_. The Poincare map T|n+ must 
be symmetric to r|n_; so we have u± = 0, u+ = Rv-, B+ = RB- and A+ = — A_ = A in 
equation (l5|) (the expanding direction u is aligned with the vector ei). 

Separatrix value. Assume that the separatrix value A satisfies the condition 0 < |A| < 2. 


The first two conditions correspond to a codimension-2 bifurcation in the class of R- 
symmetric systems. A generic unfolding of this bifurcation is given by a two-parameter 
family of i?-symmetric systems which depends on small parameters (/r, e) where p splits 
the homoclinic loops, and e is responsible for the change in the saddle value. We, thus, can 
assume fj, = u+ and = 1 — e in ([5]), and the rest of the coefficients will be continuous 
functions of {p,e). Thus, the Poincare map T will have the form 


u = f{u,v) = {p + A\u\^ +p{\u\,Qs{u)V,l^,£))s{u), 

V = g{u,v) = Qs{u)iv++ + ql\u\,Qs{u)V,P,£)), ^ 

where s{u) is the sign of u, and Qg denotes identity if s = 1 and the involution R (restricted 
to the u-space) if s = —1. The functions p,q satisfy 


p,q, 


d{p, q) 

dv 


o{\u 




d{p, q) 

du 


o{\u\ ^); 


( 8 ) 


see [36], Section 13.8. 

The following theorem is a particular case of the criteria for the birth of the Lorenz 
attractor which were proposed in mn- 


Theorem 2 In the plane of parameters {p, s) there is an adjoining to zero region Vla such 
that the system possesses a Lorenz attractor for all {p,e) G Vla- 


Proof. It is known m that \iv > 1 , then stable periodic orbits are born from homoclinic 
loops. All orbits in the Lorenz attractor are unstable, so we will focus on the region z/ < 1, 
i.e. e > 0. In order to prove the theorem, we show that for all small e there exists an interval 
of values of p for which the Poincare map ([7]) satisfies conditions ([ 6 ]) on some absorbing 
domain D. 

First of all, we check the conditions under which the invariant domain D exists. According 
to [ 6 ], the domain D lies within the region bounded by the stable manifolds of saddle periodic 
orbits that are born from the homoclinic loops at Ap < 0. Therefore, we will consider the 
case p < 0 if A > 0 and ^ > 0 if A < 0. We define 

D = {\u\ < \p\ + 6{p)} (9) 
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where 6 is some function of [i of order o(/i). Then, the forward invariance condition TD C. 
int(D) reads, according to iffl) . itHll . as 

+ 0(1^1 -^)-|/r|<M + o(M), 
or 

Ia‘I>(^(i+ ■■■))■ (10) 

where the dots stand for the terms tending to zero as ^ 0. Note that if 1^41 > 2, the 
value at the right hand side of (jlOp is bounded away from zero at small e > 0, i.e., there 
are no small ^ satisfying (|inh . Thus, we consider the case |^| < 2, so for all small e this 
inequality has solutions in any neighbourhood of the origin in the (/r, e) plane. 

Now, let us check conditions ([6]) in D. By ([7]), dl]), we have 

/n = ^(l-^)l^r^s(u) + o(|nr^), {f,gy^ = o{\u\^-^), g'^ = 0{\u\-^). 

Thus, 

\\9'uifL)-X = 0ii), ||/;i|„=o(|H'-^). 

It is immediately seen from these formulas that conditions ([6]) are fulfilled everywhere in D 
if 

H< (1^1(1+ ...))i (11) 

and /i and e > 0 are small enough. 

Conditions (jlOp . (jUP define a non-empty open region in the parameter plane near (/i, e) = 
0. By construction, the intersection of this region with {Ag, <0, e > 0} corresponds to the 
existence of the Lorenz attractior. □ 

Different versions of this theorem were proved by Robinson in [ini US] under an addi¬ 
tional assumption (an open condition on the eigenvalues of linearisation matrix at O). This 
assumption made the results of [iniiis] inapplicable to the extended Lorenz model in the 
domain where we consider it (the case of small /3 in ([3|)). We therefore gave here a proof of 
the Shilnikov criterion without the extra assumptions. 

The difference with the Robinson’s proof is that we do a direct verfication of Afraimovich- 
Bykov-Shilnikov conditions for the Poincare map T, while Robinson performs a factorisation 
of the flow by the strong-stable invariant foliation. As we mentioned, the Poincare map T 
has a C^-smooth invariant foliation The smoothness of is not important for the 

analysis of the dynamics in the Lorenz attractor from the point of view of the topological 
equivalence. On the other hand, it can be useful in the study of some statistical properties, 
like correlationd decay in the flow, etc. [38]. The smoothness is, essentially, a consequence of 
two properties, that the map T is hyperbolic and the expanding direction is one-dimensional, 
see [52]. Despite the hyperbolicity of the Poincare map T in D, the suspension flow in D 
is not hyperbolic, as the set D contains the equilibrium state O. The flow is only pseudo- 
hyperbolic [laiSS], i.e. it has a codimension-2 strong-stable invariant foliation such that 
the flow is contracting along the leaves of the foliation and area-expanding transversely to it. 
The foliation is generated by in the sense that the intersection of the orbit of any 
leaf of with the cross-section 11 consists of the leaves of One infers from this that 
the foliation is C^-smooth in 'D\0. However, in order the foliation to be smooth at 
O too, an additional open condition on the eigenvalues of the linearisation matrix at O must 
be fulfilled. This was the condition imposed in [n3S2] and the smoothness of the foliation 
was used in the Robinson proof of the birth of the Lorenz attractor in an essential way. 
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In what follows we apply Theorem[2l in order to prove the presence of the Lorenz attractor 
in system (j4]). It suffices to find parameters (a, A,/3), corresponding to the existence of a 
homoclinic butterfly with 7 = 1 , compute the separatrix value A on the homoclinig orbit, 
and check that 0 < A < 2. 


3 Lorenz attractor in system (Hj) at small (3 


3.1 Homoclinic butterfly 

At /3 = 0 system (jl]) has an invariant plane Z = 0. If, in addition, A = 0, then the restriction 
of dH to this plane is a Hamiltonian system of the form X = X — X^. At zero energy level 
it posesses a pair of homoclinic solutions x{t) = zbxo(t), yo{t) = ±xo{t), where 


xo{t) 


V2 

cosht 


( 12 ) 


Note that by setting a = 1 we make the saddle value a at the origin vanish. 

Next, near the point (a,A,/?) = (1,0,0) in the space of parameters we will find the 
bifurcation surface which corresponds to the existence of a homoclinic butterfly in system 
dH). We do it by introducing a small parameter /x = (a — 1, A, /3) and expanding formally the 
homoclinic solution in powers of y: 


X{t) = xo{t) + xi(t) + 0{\\yf), Z{t) = z,{t) + 0{\\yf) (13) 


where xi,zi = 0(||/u||). One obtains a linear system (with time-dependent coefficients) for 
the order-11/X11 corrections xi, 2 ;i; this system has a solution which tends to zero as t ^ ±oo 
when a certain linear functional of y vanishes. We check that this functional is not identically 
zero, so it vanishes on a certain plane in the /x-space. It is well known (see [5l] for a proof) 
that this implies the existence of a surface Mhom (tangent to this plane at /x = 0 ) such that 
system (jH) has a homoclinic solution of form (|13p when the parameters belong to Mhorn- 
We consider a more general system: 

X + V'{x) = f{x, x)z + g{x, x)y -I- z‘^(j)i{x, x, z, /x), 

z =-(a(fj.) + p(x,x))z + q(x,x)fi + z^^ 2 (x,x,z, fi), ^ 

where /x is an n-dimensional vector of small parameters, and V, /, g, p, q and 4>i^2 are 
certain sufficiently smooth functions. We assume that H(0) = H'(0) = 0, IA"(0) < 0, 
g{0, 0) = ( 7 ( 0 ,0) = 0, p{0, 0) = 0, and a(0) = ao > 0. By construction, the origin O is an 
equilibrium of (fT4l) for all /x. System dlj is a particular case of (fT4)l . where one takes 

^(2;) = -y + y> f = -x, P = 0, g = {0,-x,0), q = {0,0,x‘^). (15) 

When /U = 0, system (I14p has an invariant plane z = 0, and the restriction onto this plane is 
Hamiltonian: 

X + V'{x) = 0, (16) 

with the hrst integral 

x'^ -|- 2V(x) = const. (17) 

Let there exist at least one value of x 7 ^ 0 where V(x) = 0, and let x = x* he the closest 
such point to O (without loss of generality we assume x* > 0). Thus, at the zero level of 
integral (flTl) there exists a separatrix loop Lq to O lying at z = 0. We take a parameterisation 



xo{t) of To such that xo(0) = x* and, therefore, xo(0) = 0- It is clear that xo{t) tends to zero 
exponentially (with the rate uj = ^"(O)) as t —)• ±oo. 

For fi ^ 0, the plane z = 0 is no longer invariant but the homoclinic loop F^ may still 
exist if a certain codimension-1 condition is fulfilled. Note that if an additional {x —?■ —x) 
symmetry condition is imposed on (jl4p . then there will be two homoclinic loops to O. We 
will seek for a homoclinic solution in form (fT^ . where {xi{t), zi{t)) = 0(||/r||). One obtains 
from (I14p the following system 

xi + V''{xo{t))xi = f{t)zi + g{t)^, 
ii = -(a(^) +p{t))zi + q{t)n, 


where if,g,p,q){t) = {f,g,p,q){xo{t),xo{t)), i.e. they are explicit functions of time. 
System (jl8jl is an inhomogeneous linear system; its homogeneous part is given by: 

xi + V"{xo{t))xi = f{t)zi, 
zi = +p{t))zi. 

System (jl9h has the following three linearly independent solutions: 

where 


(19) 


Xi{t) = Xo{t), X2{t) = Xo{t) 


ds 

Wr 


Xsit) = -xo{t) 


t +CXD 

f dw 




xois)f{s)h{s)e °‘^ds, Z^{t) = e “ h{t), 


and 


- jp{s)ds 

h{t) = e 0 


( 20 ) 


Note that xo(0) = 0 so that, generally speaking, the integrals of the form f in these 

formulas do not converge at u) = 0. We, nevertheless, keep this notation in the following 
exact sense: 


t 


t 

/ 


F(w) 

XoM 


dw = lim < 
€—^~!“0 



m 

xliw) 

exg(O) 

F{w) ^ 

„, (dw + 

m 

xliw) 

exg(O) 

infinity as t 

—)• 0, : 


t > 0, 


t < 0. 


( 21 ) 


finite limit. 

It is easy to see that the asymptotic behaviour of these solutions is represented by the 
following table: 

t —oo t -|-oo 


e-uit 


^—at 


where to = F"(0). 


Xi{t) 

X2{t) 

X3{t),Z3{t) 
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Once we know the solution of the homogeneous system, we find the solution of the inho¬ 
mogeneous system (fTHI) : 

* F{w) 


xi{t) = Xo{t) 


Xq{w) 


dw ■ n, 


zi{t) = e 


h{s) 


q{s)ds ■ /i, 


( 22 ) 


where 


W 

F{w) = J xo{s) 


f{s)h{s)e~ 


h{v[ 


-q{v)dv -|- g{s) 


ds 


(23) 


and h is given by (12011 . 

For solution (j22ll to correspond to a homoclinic orbit, it must tend to the origin in both 
directions of time. It is easy to see that lim zi{t) = 0 and lim xi{t) = 0. However, as 

t—>-dioo t—>•—oo 

t —7> -boo, the function xi{t) will converge to zero only if 

F(-|-oo) ■ g = 0. 

As we mentioned, when this condition defines a codimension-1 subspace in the /i space, there 
exists a codimension-1 manifold Mhom tangent to this subspace at ^ = 0 such that system 
m has a homoclinic loop at fi G Mhom- Thus, the surface Mhom exists when 


F(-l-oo) / 0. 


(24) 


Now, let us apply the obtained result to system (j4]). By plugging (fT^ . (fTHI) into (l23]l we 
find that 


F{w) ■ g= f xois) 


fj,3Xo{s)e ® f x^(v)e'^dv + ^ 2 x 0 ( 3 ) 


ds = 


= --(1 -b tanm w)X + { tx — 4arctam e"' -b 


55 


15 


12 


--b 


4 cosh^ w 

4(tanh u) — 2) arctan e’" , Stanh^tc 4tanh^rc\ ^ 

H—- 7 ^ -b 4 tanh w H-:-^-^- ) /?. 


(25) 


coshtc 


This gives us 


F(-boo)/r = ( ^ - /3 - ^A. 


Condition (j24p is, obviously fulhlled, so we obtain the existence of the sought homoclinic 
butetrfly for parameter values belonging to a smooth Mhom of the form 

4 


/? = 


rA -b 0(A -b (o — 1) ). 


32 - 37r2 

The corresponding homoclinic solution is given by (I13p where, according to 
have 


(26) 


we 


Xl{t) = - 


t 

xo(i) / TaTA 

J x^(w) 


/55 


—-(1 -b tanm w)A + { tx ~ darctan^ e"' -b 
o V X Z 


15 


dftanhrc — 2 ) arctan e"' , 3tanh^t(; dtanh^tc , ^ 

H—^- 7 ^-bdtanhrc H-^-^-^- ) fd 


coshtc 


4 cosh^ w 
dw 


+ 


( 27 ) 


zi = jjL^e 


-t 


Z 

J xl{s 


s)e^ds = 2e * ( 2arctane* — 


1 


cosh t 


/?• 
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Since the condition cr = 0 can be written as A = 2(a — 1) + 0((a — 1)^), we obtain from 
(|26h the following parametrisation for the bifurcation curve corresponding to the existence 
of a pair of homoclinic loops to a saddle with zero saddle value in a small neighbourhood of 

fi = 0: 

O' = 1 — £ O(e^), A = 2e + O(e^), /3 = ^— ^£ + O(e^), e > 0. (28) 

We need to check that all conditions of the Shilnikov criterion are satisfied by the homo¬ 
clinic loops that exist when the parameters belong to this curve. First, let us check that the 
homoclinic loops at e > 0 tend to the equilibrium along the leading direction. 

Make the following change of variables in system Q in order to align the coordinate axes 
with the eigendirections of the linear part of the system at O: 

X = —au + V, Y = u + av. 


The system takes the form 


1 


u = - u + ^ „ 

a I + 


" uZ--^vZ- ^ —^3 


1 + a^ 


1 + a- 


:{v — auY 


V = av + 


a 


l + a^ 


tuZ — 


a 


l + a^ 


,vZ- 


a 


1 + a‘ 


■{v — auY 


Z = —aZ + Y{v — auY- 

Here the coordinates u and Z correspond to the stable directions, and v corresponds to the 
unstable direction. Note that u corresponds to the strong-stable direction (i.e. the corre¬ 
sponding eigenvalue is farther from the imaginary axis than the eigenvalue that corresponds 
to the Z direction; recall that a < 1). Thus, the equations of the local strong stable and 
unstable manifolds near the saddle are 

W,-: Z = 7iu2 + 0(|u|3), v = 0{u^), 

WYo,-. Z = 72u2 + 0(|u|3), u = 0{v^), 


where the coefficients 71 and 72 are found by equating the coefficients of the power series 
expansion of the conditions of invariance of and, respectively, One obtains 


71 = 


f3a 

a ^-2 


< 0 , 


72 = > 0. 

6 a 


This implies that the manifold WiYY\0 lies in the region Z < 0, while lies in the 
region Z > 0. Finally, since Z > 0 everywhere at Z = 0 except for the equilibrium point, 
the global unstable manifold never crosses Z = 0, i.e. the unstable separatrices of the 
saddle can never enter the region Z < 0, hence they cannot enter Wf^^. This proves that 
both homoclinic loops enter the saddle along the leading direction - the Z-axis; as they come 
from the same side Z > 0, they are tangent to each other at t = -|-oo and form a homoclinic 
butterfly, as required by the Shilnikov criterion. 


3.2 The separatrix value 

To verify the last condition of Theorem [2] we will compute the separatrix value A for small 
values of /3 and check that it lies in the range 0 < |H| < 2 which will finalize the proof of 
Theorem [TJ In order to determine the separatrix value we will use the definition from |34j^ . 

^Note that different but equivalent definitions can be found, for example, in MISS] 
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Let the equation of the homoclinic loop be u{t) = {x{t),y{t), z{t)), t G (—oo,+oo). We 
consider the linearisation of ([4]) near this solution: 


? = cm 


(29) 


Cit) = 


where 

0 1 0 

1 — 3x'^{t) — z{t) a — Xja —x{t) 

2/3x{t) 0 —a 

For any two vectors ^i(t) and ^ 2 {t) satisfying (f2^ their vector product rj = ^ 1 X ^2 evolves 
by the rule 

7 ) = D{t)r] (30) 


with 


D{t) =trC(t) •Id-C^(t) = 


As x{t) and z{t) tend to zero when t 
is determined by the limit matrix 

= 


— 1/a 3x^{t) + z{t) — 1 —2/3x{t) 

— 1 —a 0 

0 x{t) a — 1/a 

-)• ± 00 , the asymptotic behaviour of solutions of ((30]) 

-1/a -1 0 

-1 -a 0 

0 0 a — 1/a 


The eigenvalues of D^o are 0, — (a + 1/a) < 0 and a — 1/a < 0. Thus for t —+00 each 
solution of (l30t) tends to a constant times the eigenvector y* = (—a, 1,0)^ corresponding to 
the zero eigenvalue. At the backward time t —>■ —00 all solutions grow to infinity except for a 
one-parameter family of solutions which tend to rj* multiplied to some constant. Thus, there 
is only one solution 7/(t) such that lim r]{t) = r]*. For this solution, we have 

t^ — OO 

lim r](t) = An*. 

t^+00 

The constant A here is indeed the separatrix value; its absolute value is the coefficient of 
contraction/expansion of infinitesimal areas near the homoclinic loop, and its sign indicates 
the orientation of the loop. It is clear that 


|A| = sup lim 

t^+00 


\vi-t)r 


(31) 


where the supremum is taken over all the solutions of (|30p . 

To compute the separatrix value, perform a coordinates rotation such that the eigendi- 
rections of the matrix become the coordinate axes. The system (I3UI) takes the form: 


a 


1 + a- 


vi = - 


V2 = - \a + 


:{3x^{t) + z(t))(vi + V 2 ) + 


2af3x{t) 
1 + a^ 


-V 3 


a 


V2 + 


a 


1 + a- 


■(3x^(f) -h2(t))(ui +V 2 ) - 


2af3x{t) 
1 -(- a^ 


-V 3 


(32) 


7)3 = ^a - — J V3 + x{t)vi + x{t)v2. 

Here vi is the coordinate in the direction of the eigenvector r]*, corresponding to the zero 
eigenvalue of the limit matrix. Thus we are interested in the solution of ((32 p satisfying the 
conditions 7;i(— 00 ) = 1, 7;2(— 00 ) = U 3 (— 00 ) = 0, and the separatrix value is A = 7;i(+oo). 
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Introduce a new variable Vi = v\ + V 2 to simplify the system: 


V4 = - 



V 2 


V2 = - 



1 + 


2af3x(t) 

~r ,—^^3 
1 + 


(33) 


V3 



V 3 + x(t)v 4 ; 


the solution we are looking for satisfies n 4 (—oo) = 1 , V 2 {—oo) = ?; 3 (—oo) = 0 . As z; 3 (+oo) = 0 
automatically, we have that A = U 4 (+oo). 

It is not hard to show that system (l33]l has a unique solution satisfying these conditions 
for all /3 and a G (0,1]. Just note that the problem is equivalent to the following system of 
integral equations: 


V4{t) = 1 — i a + 


a 


V2{s)ds, 


V2{t) = j 




a 


1 + a- 


r(3x^(s) + z(s))u4(s) 



e(a *^x(s)u4(s)ds, 


— OO 


2afix{s) 

\ + 




ds, 


(34) 


and, since x(s) and 2 ;(s) tend exponentially to zero as s —>■ —oo, it immediately follows that 
the integral operator on the right-hand side of this system for t close to —oo is uniformly 
contracting in the space of bounded continuous functions on (—oo,t] with appropriately 
chosen exponential weights for V 2 and V 3 . Once the solution is shown to exist up to a certain 
value of t, it is continued to all larger values of t as a solution of the Cauchy problem for 
system (I33p . Since the solution is obtained as a of a contracting linear operator, which 
depends smoothly on parameters, the solution also depends smoothly on the parameters. 
Because of the uniform convergence of the integrals in (I34p as t ^ -|-oo, the separatrix 
value A = U 4 (-|-oo) is also a smooth function of the parameters, i.e. it can be found by an 
asymptotic expansion. 

Thus, we will seek for the solution as a power series in e using formulas (I28p . The unknown 
functions are represented as 

Vi = Vi + vjs + i = 2,3,4. 


The equation of the homoclinic loop has the form 

x{t) = xo{t) + Xi{t) + O(e^), y{t) = xo{t) + xi{t) + O(e^), z{t) = zi{t) + 0(£^). 

For e = 0, system (l33]l is rewritten as: 

V 4 = - 2 v 2 , vl = - 2 v 2 + ^xl{t)vl, V 3 = Xo(t)v 4 . (35) 

The first two equations here do not depend on U 3 and are reduced to the equation 

V4 + 2x4 ^xl{t)vl = 0. 
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Its solution satisfying ti4(—oo) = 1 is 



sinh t 
2 cosh^ t 


(36) 


from which we obtain 


vl{t) = 


1 


4 cosh^ t 


(1 — 2tanhf), 


t 


U3 (t) = —— (tt + 4 arctan tanh —h 


4v^ 


2 cosh t 


(2 — tanht)). 


Note that v^(t) vanishes as t —>■ +00. Therefore, the asymptotic expansion for the sepa- 
ratrix value has the form A = + O(e^). To complete the theorem, we need to compute 

Ai and show that Ai > 0. 

The hrst order in e terms satisfy the following system: 


vl = -2vl, vl =-2vl +-xl{t)vl +fi{t), vl = xo{t)vl +f 2 {t), 


hit) = -^iQxoit)Mt) + zi{t))vl{t) - 

f 2 {t) = 2v^{t) + xi{t)v2. 


xo(t)v^(t), 


(37) 


The first two equations do not depend on U3, so we will solve the following equation to find 


vr. 


vl + 2vl + -xl{t)vl + 2/1 (f) = 0 


(38) 


with the boundary conditions U4(—00) = 0 and U4(+oo) = Ai. Two independent solutions 
of the homogeneous part of (f38]l are 


, , sinh t f , , cosm t + 3t tanh t — 3 + 

yht) = -TTre > y 2 it) = - 


yiit) y 2 it) 


= e 


„-2t 


cosh^f ’ 2coshf 

with the Wronskian equal to 

yi{t) mit) 

so that the general solution of (I38p is written as: 

vl{t) = Ci{t)yi{t) + C 2 it)y 2 {t), 

and the coefficients C'i(f) and (72 (f) are determined by the following formulas: 


t t 


C'i(t) = 2 / y 2 {s)fi{s)e^^ds, C 2 {t) =-2 I yi(s)/i(s)e^*(is. 


(39) 


(40) 


It remains only to calculate the limit Ai = lim {Ci{t)yi(t) + C 2 it)y 2 it)). We have 

t—>-+oo 


lim Ci{t)yi{t) 

t —^“)~00 


f 2/2(s)/i(s)e^®ds 


lim 

t —^H~oo 


—oo 


1 

yiW 


lim 

t—>-+oo 


y2(t)fi(t)e^* 

nit) 


= 0. 
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H-oo 

As lim y 2 {t) = —1/4, this gives us Ai = 1/2 f yi(t)fi{t)e'^^dt. To compute this integral, 

-oo 

we split it into three parts. Taking into account that yi{t)xQ(t)v\{t)e^^ = -{2x‘f^{t) — a:Q(t)) 

8 

we obtain 


+ 0O +00 t 

J yi{t)xo{t)xi{t)vi{t)e‘^^dt = ^ j {2xl{t) - xl{t))xo{t) j 


3 f x^{t) r F{w) 

' xliw) 


8 

and 


dw 


t=-\-oo 


t=—oo 


3 

48 


+00 


3x^{t) - x^jt) 

xl{t) 


F{t)dt = - 


16 


9(32 - 37r2) ^ 


+ 00 




— OO 
+ 00 


— [ yi(t) - ^—^xo(t)uo=—I - 

J ^ ^^32 - 37r2 ^ J 3 32 - 37^2 

— OO 

This result gives us the following formula for the separatrix value in a small neighborhood 
of point y = 0: 


A = 


TT 


+ O(e^). 


32 - 37r2 

Thus, for all small e > 0 we have 0 < A < 2 which means that the orientable Lorenz 
attractor is born from the homoclinic butterfly. Theorem [1] is proven. □ 


4 Discrete Lorenz attractors in three-dimensional 
diffeomorphisms 

In this section we will apply Theorem [T] to prove the birth of a discrete Lorenz attractors at 
certain codimension-4 bifurcations of three-dimensional Henon-like maps. Consider a map 
(x, 7/, z) (x, y, z) of the form 

x = y, y = z, z = Bx + f{y,z,£), (41) 

where / is a smooth function, e is a set of parameters, and the Jacobian i? is a constant. 
Fixed points of the map are given by the equations 

X = y = z, x(l — B) = /(x,x). 

The characteristic equation at the fixed point x = y = z = xq is 

\^ - A\^ -C\- B = 0, 

where A = //(xo,xo), C = fy{xo,xo). When {A,B,C) = (—1,1,1) the hxed point has 
multipliers (j-l, —1, —1). We will study bifurcations of this point. To do this, we shift the 
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origin to the fixed point and introduce small parameters ei = l — B, 82 = 1 — C, and 
= —1 — A. The map will take the form 


x = y, y = z, z = {1 - ei)x + {1 - e 2 )y - {1 + e 3 )z + ay'^ + byz + cz'^+ 
+diy^ + d2y‘^z + + <^4^^ + ..., 


(42) 


where the dots stand for the rest of the Taylor expansion. 

In [12] there was shown that if the following inequality is fulfilled 

'll) = {c — a){a — b + c) > 0 , 


(43) 


then map (|42l) near the fixed point at zero satisfies conditions proposed in |3|, which implies 
that the second iteration of the map is close, in appropriately chosen rescaled coordinates, 
to the time-1 shift by the orbits of the Morioka-Shimizu system. In other words, the map 
is a square root of the Poincare map for a certain small, time-periodic perturbation of the 
Morioka-Shimizu system. As the Morioka-Shimizu system has the Lorenz attractor, its small 
time-periodic perturbation has a pseudohyperbolic attractor [53] which corresponds to the 
discrete Lorenz attractor in the root of the Poincare map. Thus, we obtain that a discrete 
Lorenz attractor containing the fixed point exists in a small neighborhood of zero for a certain 
region of small parameters £1,2,3 in map (I41|l . 

One can easily repeat the calculations from m for 1/^ < 0 and obtain the same normal 
form in this case, but with one difference - the scaling factor for the time will be negative, 
which means that the presence of the Lorenz attractor in the Shimizu-Morioka model implies 
the presence of a discrete Lorenz repeller in map (I4ip in the case i/j < 0. We mentioned that 
the existence of the Lorenz attractor in the Morioka-Shimizu model is definitely true but not 
formally proven, which means that we are confronted with the same problem when trying to 
establish the existence of the discrete Lorenz attractor or repeller in the 3D Henon map. 

We bypass the problem by considering the case of an additional degeneracy; namely, we 
study bifurcations when -0 vanishes. The list of rescaled normal forms corresponding to a 
hierarchy of a certain class of degeneracies of the (-1-1, —1, — l)-bifurcation is obtained in [5]. 
It was shown there that the extended Lorenz model does appear as a normal form in some of 
the degenerate cases. Therefore we can apply Theorem 1 and, thus, obtain an analytic proof 
of the existence of a discrete Lorenz attractor in map (142 p for some open region of parameter 
values; see Theorem [3| below. 

There are two multipliers in the formula for so we consider two possible case£|: 


Case I : (a — 6 -I- c) = 0, (o — c) 0 
Case II : (a — 6 -I- c) 0, (a — c) = 0. 


(44) 


We introduce the fourth independent small parameter £4 as follows: 


Case I : £4 = -(a — 6 -|- c) 
^ 1 

Case II : £4 = -(a — c). 


(45) 


Consider a four-parameter family of maps of type (I42p which depends smoothly on £^ 2 3 4- 
In particular, the coefficients a, b, c, and di depend on £^ 2 3 4 smoothly. The following lemma 
provides the normal form for this codimension-four bifurcation in both cases: 

"‘Note that the situation when both {a — b + c) =0 and (a — c) = 0 has a higher degeneracy, i.e. we would have 
to deal with codimension at least five here. We do not consider this case in this paper. 
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Lemma 1 For any sufficiently small £12 3 a second iteration of map m is close, in 
appropriately chosen rescaled coordinates, to the time -1 shift by the orbits of: 

X = Y, 

Y = X - \Y - XZ ±X^ + ... (46) 

Z = -aZ + /3X‘^ + ..., 

in Case I, and 

X = Y, 

y = X - Ay + / 3 XZ + RiX^ + R 2 YZ + RsXZ'^ + ... (47) 

Z = -aZ + x2 + z2 + .... 

in Case II. 

Here the dots stand for vanishing at s = 0 terms, the parameters a, ft, and A are functions 
of e which can take arbitrary finite values, and the coefficients Ri{e) are bounded as small 
parameters £ vary. 

As we will show below, the sign of the X^ term in formula (1461) coincides with the sign of 

G = 4(di -d 2 + d 3 - di) - (5 - 2cf. (48) 

Therefore, when G < 0, the normal form (I46|) for Case I coincides with system ([1|, which 
possesses a Lorenz attractor near (a,/3,A) = (1,0,0) according to Theorem [H This fact, by 
Lemma m, implies the existence of a discrete Lorenz attractor in the original map (I41jl . 
Namely, the following theorem is valid: 


Theorem 3 Assume that for some £ = Eq map possesses a fixed point M{xo,xo,xo) 
with multipliers (1,—1,—1), such that 


a — 6 + c = 0, G<0, 


(49) 


where G is given by formula and a,b,c, di are the coefficients in the Taylor expansion 
Then in some small neighborhood of Eq in the parameter space there exists a domain 
Vdla such that map has a discrete Lorenz attractor when £ G Vdla- 

It remains only to prove Lemma [H 

Proof. Let iz be the multiplier of the zero fixed point of map (1421) which is close to 1: 

El + £2 Y £3 


v = l- 


+ 0(||e|| 


Perform the following linear coordinate transformation: 

ui = -ux + y, U2 = nx + {u-l)y- z, U3 = n{l - £i)x + {1 - £i + n{l - £2))y+ n‘^z. 
Map ()42p takes the following form 


Ul = -Ul - U2, 

U 2 = -hiui + (-1 + 52)u 2 - —{a{2ui +U 2 + ^ 3 )^ - h{2ui +U 2 + U 3 )x 

1 

x(2tti + 3rt2 - U3) + c{2ui + 3u2 - tts)^) + -—{-di{2ui + U2 + it3)^)+ 

d4 

+^ 2 ( 2^1 + 3?X2 - U3){2ui +U2 + U 3 Y - d3{2ui + 3ti2 - U3)^(2ni + U 2 + ^ 3 )+ 
+d4{2ui + 3u2 - us)^) + 0(||u|p||e|| + ||tt||'^), (50) 

-U 3 = (1 - 63 )u 3 + —(a( 2 rii + tt 2 + ^ 3 )^ - b{2ui + U 2 + U3){2ui + 3u2 - 1 ( 3 )+ 

1 

+c{ 2 ui + 3 U2 - U3f) - -—{-di{ 2 ui +U2 + Us)^) + d2(2ui + 3U2 - U3)x 

o4 

X (2tti + tt2 + U3)'^ - ^3(2^1 + 3 U 2 - tt3)^(2tti +U2+ U3) + 

+d4^{2ui + 3u2 - us)^) + 0(||u|p||e|| + ||ti||^). 
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where 


(5i = + 63 


1 — 61 — £2 + £3 


(52 = 1 - zy - 63 = 


V 2 

61 + 62 — 863 


+ o(lkf), 


+ 0(||6|| 


X 1 £1 + £2 + £3 , ^/|| || 2 \ 

(53 = 1 - Z/ = - - - + 0(||6|| ). 

Note that the linear part of (|50p is in the Jor(dan form at 6 = 0. 

With a close to identity polynomial change of coordinates we kill all non-resonant quadratic 
and cubic terms. One can see that the system takes the following form after that: 

Ul = —Ul — U2, 

1 1 

U 2 = -(5irti -I- (-1 -I- S 2 )u 2 - -(a - c)uiU 3 - -(a - b - 3c)u2U3 + B 3 oouf + Bio 2 Uiu§-l- 

-09(||rz|p||6|| ||n|p|tZ 2 | + ||u|n, 

U 3 = (1 — <53)113-1 -—(a — 6-|- c)ui -j- —— (fl-|-6 + c)u3 -|- 0(||ii||^||6|| -|- |iiiit2| ~l~ '^2 
4 ib 

(51) 


where 


B300 — + b'^ + 5c^ + ab — 56c) -|- — (—di + <^2 — ds -|- d^), 

1 1 ° 

Hi02 — ^ ( 2 a^ — b"^ db ‘Zac — 6c) -t- 737(— 3 di — d2 d3 8(^4). 


>102 V—— - , ,__ , , 

iU2 ^ 28 ' ^ 32 ' 

Next, one checks that the second iteration of the map (15ip coincides, up to terms of order 
l^(ll'*^lPlkll + ll^ll^)) with the time-1 shift by the flow of the form 

ill = Pi 111 + (2 -I- P2)u2 “ t(® “ c)miU 3 -I- 7^ (a -I- 86 -I- 5c)iZ2ii3 + 0 (||n|p), 

1 1 _ _ 

U 2 = Pi 111 - P2U2 + x(o - <3)111113 - -rib + 2 c)u 2 U 3 + Hsooii? + Hio 2 itiiii + 0 (|iZ 2 | ||ii|P), 

1 1 

II3 = -P 3 U 3 + ^(<1 - ^ + <3)111 + o (« + b + c)u\ + 0(|lllll2| + ni + \\u\f), 


where B 300 = G/128, 




(52) 


and 


so 


Pi 1 + P2 
pi -P2 


= In 


1 1 
61 1 — 62 


Pi = 


£1 - £2 + £3 


+ 0(||£|| 


P2 = 


7 £i - £2 - 563 
24 


+ 0(||6|| 


Pi = 


£1 — £2 + £3 


+ 0(||£| 


P2 = 


£1 - £3 


+ o(lkl 


2 - 2 

Introduce new parameters and time via the following formulas: 

t = t/s, s = s/^, a = p 3 /s, A = P 2 /S. 

Now we consider Cases I and II separately. 

In Case I with 64 = — (a — 6 -|- c) we perform the following scaling of the coordinates and 
parameter 64: 

Ul = i^iX, U 2 = 122 Y, U 3 = Z/3Z, /3 = -4>f(a - c) —, (58) 

s 
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where 


(54) 


1^2 ! — 2 ^ 

X > 0, Ui = [2 + P 2 ) — , V2 = =-• 

s a — c 


Then, equation ([52]1 recasts as 

X = Y + Ois), 


kG 


Y = X-XY-XZ + —x^ + 0(s) 

16 ^ ^ 
Z = -aZ + /3X2 + 0(s). 


(55) 


According to [5], the equality G = 0 is the higher order degeneracy condition, which we 
do not consider here, i.e. we assume that G 7^ 0 when e = 0. After setting x = 16/G, system 
(f55|) takes the form ()46)) up to 0(s)-terms (recall that x > 0, which does not allow us to 
control the sign of the A^-term). 

In Case II we have S 4 = —{a — c) and the scaling of the coordinates and the fourth 
parameter is performed as follows: 


Ml = viX, U2 = V2Y, M3 = M3Z, 13 = 2x(a — b + c) — 

where 

Ml = (2 + P2) —, M2 = •v/x'S^, z ^3 = -b + c)s. 
s 

In this case system (I52p takes the form 

X = Y + Ois), 

xG 1 

Y = X-XY + ( 3 XZ + -- -x(6 + 2c)(a - 6 + c)YZ+ 

IQ 4 V 

+2x^(a — 6 + c)^ B1Q2X Z^ + 0{s + £4) 

Z = -aZ + A^ + -b + c)^Z^ + 0 (s). 


(56) 

(57) 


(58) 


Finally, taking x = --;-^ we obtain formula (147)1 . Lemma is proven. □ 

(a — 0 + c)^ 

Now, according to Lemma (H the flow normal form of map (|41ll for the codimension-four 
bifurcation under consideration coincides with system ([4]) up to asymptotically small terms. 
This system has a Lorenz attractor in some domain of the parameter space due to Theorem [H 
and this implies the existence of a discrete Lorenz attractor in map (1411) . Theorem[3]is proven. 
□ 

We remark that system ()47)1 was not studied before, so the question of the existence 
of Lorenz or other attractors in this system is open. However, a similar system of ODEs 
was obtained in |55j as a finite-dimensional reduction of the Gray-Scott PDE model with a 
codimension-two singularity. Numerical experiments with the obtained ODE model revealed 
a Lorenz-like chaotic behaviour |56] . Thus, the analytic study of Lorenz attractors in the 
normal form ()47p can provide results relative to the spatio-temporal chaos in heterogeneous 
media. 
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